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Abstract. We show that the third Goldstone mode in the two-species 
condensate mixtures, which emerges at phase-separation, gets hard¬ 
ened when the confining potentials have separated trap centers. The 
sandwich type condensate density profiles, in this case, acquire a side- 
by-side density profile configuration. We use Hartree-Fock-Bogoliubov 
theory with Popov approximation to examine the mode evolution and 
density profiles for these phase transitions at T = 0. 


1 Introduction 

Two-species Bose-Einstein condensates (TBECs) are uniquely different from the single¬ 
species BEG as they can be either in miscible or immiscible phase depending upon 
the interaction strengths [T]. It is, however, possible to tune a TBEC from miscible 
to immiscible domain or vice-versa in experiments [5] through Eeshbach resonance. 

In the last decade or so, TBEGs have been realized experimentally in mixtures of two 
different alkali atoms [211], or two different isotopes [5] and hyperfine states |5]. 
Such an unprecedented control over experiments has motivated numerous theoreti¬ 
cal investigations on the stationary states m, dynamical instabilities and 

collective excitations mm of TBEGs. 

The present work describes the development of a gapless m Hartree-Eock-Bogoliubov 
theory with Popov (HEB-Popov) approximation [14] for TBECs. We apply it to ex¬ 
amine the spectra of the quasi-particle excitations of trapped TBEC in separated and 
non-separated trap centers. In the first part of the work we study the hardening of 
the Goldstone mode with the increase in the separation of the trap centers in which 
the TBECs have been confined. This has profound experimental implications since, 
in experiments, the trap centers never coincide. This is due to gravitational sagging, 
and deviations of the trapping potentials from perfect alignment. In the later part 
of our study, we show the softening of Goldstone mode with the variation in inter¬ 
species interaction strengths, albeit, with co-incident trap centers. Previous studies 
have shown the existence of an additional Goldstone mode in the excitation spectrum 
at phase-separation for the transition from miscible to symmetry-broken or side-by- 
side density profiles of TBEC |1211511b] . We, here, show that in the the changeover 
from miscible to sandwich type density profile, where one of the species is enveloped 
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by the other, is accompanied by the appearance of a third Goldstone mode at the 
transition point. Even at higher inter-species interactions, the mode continues to be 
a Goldstone mode. 


2 Theory of two component BEC 

To provide a coherent description of the HFB-Popov theory for a quasi-lD TBEG, 
we start with the theory for a single-species BEC. In a cigar-shaped harmonic trap 
for which the radial frequency u!±{u!x = = w_l) ^ the axial frequency the 

time-independent ID generalized GP equation with the HFB-Popov approximation 
is given by m 

h(j) + U [ric + 2fi] (j) = 0, (1) 

where h = {—ti^/2m)d'^/dz^ + V(z) — fj,. The strength of the repulsive contact in¬ 
teraction is represented hy U = {aX)/m with a as the s-wave scattering length, and 
A = S> 1 quantihes the anisotropy of the trapping potential. The atomic 

mass of the species is given by to, and /i is the chemical potential. The quantities 
nc{z) = |(?!)(z)p, n{z) = {z,t)tp{z,t)) , and n{z) = ndz) -\- n{z) in the above equa¬ 

tion represent the condensate, non-condensate, and total density, respectively along 
the axial direction of the trap. In the equation (j){z) is the order parameter and '4i{z) 
is the field operator associated with the fluctuation. The non-condensate or ther¬ 
mal component n{z) is computed from the solutions of the Bogoliubov-de-Gennes 
equations 

{h + 2Un)uj — Ucpi^Vj = EjUj, 

-{h + 2Un)vj + Ucp^'^Uj = E^Vj. (2) 

Here Uj’s and Uj’s are the Bogoliubov quasi-particle amplitudes, and EjS are the 
eigen energy of the jth quasi-particle excitation. Applying the same formalism, the 
ground state of a quasi-lD TBEG is described by a pair of coupled generalized GP 
equations 

hkfpk + Ukk {rick + 2hk) 4>k + Ui2n^-k4>k = 0. (3) 

Here hk = {—h‘^l2mk)d‘^fdz^ + Vk{z) — fXk and k = 1,2 is the species index. The 
quantities Ukk = {akkX)/mk and 17i2 = (ai 2 A)/( 2 TOi 2 ), are the intra- and inter¬ 
species interactions, respectively, with toi 2 = toiTO 2 /(toi -I-TO 2 ) as the reduced mass. 
For the present study, we consider repulsive interactions, that is akk,<ii 2 > 0. The 
(/)fcS are the order parameters with nck{z) = \4>k{z)\'^, fik{z) = t)'0fe(z,t)), and 

nk{z) = nck{z) + hk{z) as the local condensate, non-condensate, and total density, 
respectively of the fcth species. Like in single species case, let '4’k{z)s represent the 
fluctuations, then, the non-condensate density hfc(z) is obtained from the solutions 
of the coupled Bogoliubov-de-Gennes equations 

ClUij — UiiCpIvij + [/i2(/'l {4>2'^2j — 4>2V2j) = EjUij, 

EiVij + — 17i2^i {4>2V2j — 4‘2'^2j) = EjVij, 

E2U2j — U22(p2ri2j + Ui24>2 — (plVij) = EjU2j, 

1^2X2] + U22(t>2^U2j - Ui2(t>2 ('/'iDj “ (plUlj) = EjV2j, (4) 

where = {hi + 2Uiini + Ui 2 n 2 ), ^2 = {h 2 + 21722^2 + Ui 2 ni) and tj. = —tk- 
Here m^s and VkS are the Bogoliubov quasi-particle amplitudes, and j is the energy 
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Fig. 1. The evolution of the low-lying quasiparticle eigen frequencies in the Rb-Na TBEC; 

= 10^. (a) Shows the evolution of the low-lying quasiparticle excitations as a 
function of zo, trap center separation, in the domain 0 ^ zo ^ 3.8aosc(Rb) for = lOOao, 

(b) Shows the evolution of the low-lying quasiparticle excitations as a function of flNaRb in 
the domain 0 ^ ONaRb ^ 150ao for zo = 0. 


eigenvalue index. Taking the quantum fluctuations into account, the order parameters 
c/>kS are calculated self-consistently using Eqns. mil] as elaborated in Ref. m- In the 
present study, we examine the ground state of the TBEC and compute the energy 
eigenmodes at zero temperature both in miscible (17i2 < VU 11 U 22 ) and immiscible 
{U 12 > VUiiU^) regimes with co-incident and non co-incident trap centers. 


3 Results and Discussions 

3.1 Mode hardening in ®^Rb-^^ISIa TBEC at T = 0 due to displaced trap centers 

In experiments, the gravitational potential of Earth and tilts in the external trapping 
potentials displace the minima of the traps. The potentials are no longer co-incident 
and can be replaced by the following effective potentials 

I4(z, Zo) = [z -b (-l)''z:o] ^ , (5) 

where 2zo is the separation between two trap centers, and the other symbols have their 
usual meaning. It is, therefore, pertinent to examine the variation of energies of the 
quasi-particle excitations and the topological deformation of the density profiles with 
the change in the separation of the trap centers. For our present study, we consider a 
quasi-lD Rb-Na mixture in the immiscible regime with W 2 (Rb) = 27rx4.55Hz, ^^(Na) = 
27r X 3.89Hz and w_L(Rb) = 27r x 40.2Hz, a;_L(Na) = 27r x 32.2Hz and A^Na = .^Rb = lO"*- 
Let Rb and Na be the first and second species, respectively with an = ORbRU = lOOoo, 
022 = ONaNa = 50ao, and ai 2 = ONaRb = lOOao where oq is the Bohr radius. When 
the trap centers are co-incident, that is, zg = 0, the spectra is characterized by three 
Goldstone modes. The condensate density profile assumes a sandwich geometry, in 
which Rb condensate is at the center and flanked by Na condensate at the edges. 

For Zg > 0, separated trap centers, with the breaking of the z-parity of the system, 
the energy of the second Goldstone mode of the Na condensate gradually increases 
and gets hardened at a critical value of zg. The Na Kohn mode which is transformed 
into the third Goldstone mode of the system, second for the Na condensate, at phase 
separation with z = 0, regains energy when zg > 0 to emerge as the second Kohn 
mode of the system. This is evident from the mode evolution as shown as a function 
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Fig. 2. The evolution of the low-lying quasiparticle eigen functions as a function of zq 
in the Rb-Na TBEC in the domain 0 ^ zo SS 3.8aosc{Rb) for (a) zq = 1.5aosc(Rb), (b) 
Zo = 2.5aosc(Rb)- 

of Zq in Figs. a), and the profile of the hardened mode mode is shown in Fig. 
The energy ofthe Rb Kohn mode, on the other hand, remains unchanged even in a 
separated trap setting. Furthermore, with the change from co-incident to non con- 
incident trap center, the sandwich type density profile gradually changes into a side- 
by-side density profile. The excitation spectra of the system is now identified by two 
Goldstone modes and two Kohn modes. It is to be recalled here that, throughout the 
transition, the TBECs are always phase separated. It must, however, be emphasized 
that the evolution of the Goldstone mode when the TBEG undergoes a transition from 
miscible to immiscible with zq = 0 is different from the current case. We reported 
this in our previous work m for Rb-Cs TBEC, and in the present work we discuss 
the same for the Rb-Na in the following section. 

3.2 Third Goldstone mode in ®^Rb-^^Na TBEC at T = 0 

The transition from miscible to immiscible phases in TBECs occurs when U 12 > 
\/UiiU 22 - For the Rb-Na TBEC, experimentally, it is possible to steer the system 
from miscible to immiscible domain through the Rb-Na Feshbach resonance [18]. 
This motivates us to examine the variation of eigenenergies EjS of the quasi-particle 
excitations at, and around the point of phase separation. We thus vary in the 

domain 0 ^ ONaRb ^ 150ao to study the role of inter-species interactions in the Rb-Na 
TBEC, more precisely, the effect on the excitation spectra. In the absence of inter¬ 
species interactions, that is = 0, the excitation spectrum has two Goldstone 

modes, and two Kohn modes with eigen energies and 0.85/ia;2(Rb) for Rb and 

Na species, respectively. Further more, the modes of the two species are distinct and 
uncoupled. 

For > 0, the eigen spectra of the two species are not independent, and as a 
result the modes intermix. The energy of the Rb Kohn mode remains unchanged with 
the variation in The energy of the Na Kohn mode, however, decreases with the 

increase of , and goes soft at a critical value of when the phase-separation 
occurs. This generates a new Goldstone mode of the Na BEC in the Bogoliubov 
excitation spectrum as shown in FigQb). This happens because at phase-separation, 
the miscible density profiles transform to sandwich type density profiles, with the ®^Rb 
condensate located at the center and ^^Na condensate at the edges. Thus, effectively 
the Na condensate consist of two topologically disconnected BECs and there are 
two Goldstone modes with the same |MNa| and |?^Na|, but different phases [T7]. For 
the other case, transition from miscible to side-by -side density profiles, the Kohn 
mode energy of one of the species decreases, and goes soft at phase-separation. This 
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breaks the z-parity symmetry of the system, and the zero-energy mode then starts to 
regain energy after phase-separation m- Detailed discussions on Goldstone modes 
and bifurcations in TBECs with or without soliton are discussed in our previous 
works [TTrra] . 


4 Conclusions 

The third Goldstone mode in TBECs with sandwich profile at phase separation gets 
harden when the trap centers are separated by a critical distance. This is accompanied 
by the topological change from sandwich to side-by-side condensate density profiles. 
This result has important experimental implications, and demonstrates why it is a 
major challenge to obtain sandwich type density profiles in TBEC experiments. The 
appearance of a third Goldstone mode in TBECs with sandwich type density profile at 
phase-separation in Rb-Na is another finding, which confirms and adds to the similar 
result reported in our previous work ffU- This Goldstone mode is associated with 
the fragmentation of the condensate species at the edges, Na in the present work, 
into two distinct condensates. So, effectively, the TBEC at phase separation with 
sandwich geometry, in quasi-lD, is equivalent to three coupled condensate clouds, 
and thus three Goldstone modes emerge in the excitation spectrum. 
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